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Abstract 

In this paper we analyze the relaxed form of a shape optimization 
problem with state equation 

— div (a(x)-Du) = f in D 
boundary conditions on dD. 

The new fact is that the term / is only known up to a random pertur- 
bation £(x, uS). The goal is to find an optimal coefficient a(x), fulfilling 

the usual constraints a < a < (3 and / a(x) dx < m, which minimizes 

J D 

a cost function of the form 

j (x, ui, u a (x, u))) dxdP(uj). 

n Jd 

Some numerical examples are shown in the last section, to stress the 
difference with respect to the case with no perturbation. 

1 Introduction 

The field of shape optimization problems received in the last years a partic- 
ular attention from the mathematical community, also in view of the many 
possible applications in high-tech instruments and structures, where increas- 
ing the performances or decreasing the weight, even by a small percentage, 
could be crucial. Several books on the field have been written, exploring 
the various aspects (theoretical, numerical, modelling, . . . ) that intervene 
in this very rich subject; we quote for instance pQ, [3], [1], [6], [9], [11], [15] . 
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The general framework of a shape optimization problem is the following: 
given a bounded domain D of R d and a given right-hand side /, for every 
subdomain A C D a PDE 

E A u = f 

is considered, with given boundary data. The PDE above produces a unique 
solution ua which, inserted into an integral cost function, provides the final 
cost 

F{A)= \ j(x,u A (x),Du A (x)) dx. 
Jd 

The shape optimization problem, under a volume constraint on the class of 
admissible choices, is then 

min{F(yl) : A C D, \A\ < m). 

Due to a strong instability of the class of domains, very often an optimal 
shape does not exist, and the optimization problem is usually relaxed into a 
more treatable form, where the main unknown is the coefficient of an elliptic 
PDE on the whole set D. In the present paper we consider the simplest case, 
where the PDE is of a linear elliptic type 

- div (a(x)Du) =f in D (1) 

with the boundary conditions on dD of Dirichlet type 

u = uq on dD. 

The new fact is that the right-hand side / in Q is only known up to a 
random perturbation; more precisely, if (O, 3", P) is a probability space, we 
assume that 

f(x,co) = f(x) + £(x,u), 
where the random perturbation £ is such that 

w) dP{<jj) = for a.e. x S D. 

There are few references in the literature of this sort of random or 
stochastic optimal design problems; a general T-convergence framework was 
introduced in [7], while an optimal design problem in a finite dimensional 
setting was considered in [2]. 

The homogenization method (see pQ , [19] ) and the classical tools of non- 
convex variational problems (in particular, Young measures, see [16], [17] ) 
are the two mostly used approaches in the mathematical literature to ana- 
lyze optimal design problems. We will use the Homogenization Theory in 
order to obtain the existence of a solution and some necessary conditions of 
optimality. 

In the last section we consider some simple cases of loads f(x,u) and 
we perform a numerical analysis of the optimal configurations, showing the 
differences between the deterministic case f{x) and the perturbed one f(x) + 
£(x,u). 
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2 The optimization problem 



We consider a bounded open set D C M. d with a Lipschitz boundary, two 
constants a and j3 such that < a < (3, and a given value m G f3\D\). 

We also consider a probability space (f2, 3", P) and a random map / : 
U — > L 2 (D) that we write as 

f(x,uj) = f(x) 

where £ has the property 

/ i(x,uS) dP(uj) = for a.e. x G £>. 
For every coefficient a(x) verifying 



a < a{x) < f3, 
we consider the linear elliptic PDE 



a(x) dx < m 



D 



— div (a{x)Du) = f(x, uj) 
u = on 



(2) 



which provides a unique solution u a (x,u)). Finally, we consider a cost func- 
tional of the form 



F(a) 



n J d 



j[x,uj,u a (x,uj)) dx dP(uj) 



(3) 



where j{x,uj,u) is measurable, l.s.c. in n, and such that for suitable c > 
and A G L^D x 0) 

j(x, uj, u) > A(x, — c\u\ 2 y(x,uj,u). 

More general cost functionals, of the form 



F(a) 



n J d 



j(x,uj,u a (x,Lj), Du a (x,uj)^ dx dP(io) 



could also be considered, but we limit ourselves to the simpler case 
having in mind the energy 

j(x,uj,u) = -f{x,u))u 

and the compliance 

j(x,u,u) = f(x,u)u. 
The optimization problem we consider is 



mm 



|f(o) : a < a(x) < j3, J a{x) dx < m|. 



(4) 
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Note that the optimal coefficient we look for is deterministic, that is it does 
not depend on the random variable oj. 

Besides to problem Q we will consider, mainly for the numerical pur- 
poses, the penalized one 

min |p(a) + A J a(x) dx : a < a{x) < /3 j (5) 

where A is the Lagrange multiplier of the mass constraint. 

3 The state equation 

The approach we follow in order to prove the existence of solution of the 
optimization problem @ consists in checking that our problem can be seen 
as a relaxed optimal design problem of another auxiliary optimal design 
problem and from this relaxed character we deduce the existence of a solu- 
tion. We focus on the Homogenization Method. Throughout this section, 
we denote by Xn G L°°(D; {0, 1}), n = 1, 2, a sequence of characteristic 
functions and A n € M dxd a sequence of tensors of the form: 

A n {x) = axn(x)ld + /5(1 - Xn{x))h 

with < a < (3. 

3.1 The Homogenization Method 

The homogenization method is based on the concept of if-convergence (see 
DP; 02]) [E]) [IS])- We say that a sequence of tensors {A n } n< =N -£/-converges 
to the tensor A, € L°°(D,M nx?1 ) if, for any / such that f(-,u) G H- l {D) 
P-a.e. qo € fi, the sequence {u n } of solutions of 

f -div (A n (x)Vu n (x,oj)) = f(x,oj) in D 
\ u n = on dD. 

satisfies 

J u„{-,lo) u(-,uj) in Hq(D), P-a.e. to € Q, 

\ A n Vu n (-,ui) A*Vu(-,oj) in L 2 (D) d , P-a.e. a; € ^ 

where u(-,oj) is the solution of the homogenized equation P-a.e. w G 

J — div (i„(x)Vn(x, oj)) = f(x,uj) in D, 
\ u = on dD. 

We shall write A n — > A* to indicate this kind of convergence. 
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We consider A, B (£ M dxd , {xn}neN C L°°(D,{0, 1}) a sequence of 
characteristics functions, {A n } ne ^ the sequence of matrices 

A n (x) = Xn(x)A + (1 - Xn{x))B. 

We assume (which always occurs for a subsequence) that there exist 9 £ 
L°°(D, [0, 1]) and A* G L°°(D,M dxd ) such that 

Xn (x)^9(x) inL°°(A[0,l]) 

and 

In this case A* is called the homogenized tensor obtained by the composition 
of the two phases A and B, in proportions 6 and 1 — 9 respectively, and with 
the microstructure defined by the sequence {x n }n<=N- 

In this sense the homogenized tensor A* is characterized by three com- 
ponents, the phases A and B and the proportion 9. Therefore an important 
issue is to identify all possible homogenized tensors once fixed these three 
components, this is the so-called G-closure problem. 

Fortunately, for the case of two isotropic matrices, the G-closure in the 
deterministic case is well known (see pQ, [ID], [32], |15j). We will prove that 
our "random" G-closure remains equal to the deterministic one. We denote 
Gg and Gg the G-closure associated with the deterministic and random 
equations. 

Theorem 1. Given 9 € L°°(D; [0, 1]) the G-closure of the two isotropic 
tensors aid and f3Id with proportions 9 and (1 — 9) respectively, is the set of 
symmetric matrices with eigenvalues X%, A2, A<j such that, 

X e < Xi < Xj l<i<d 

d-1 



Er 1 



< + 



. =1 Xi - a X 9 - a XJ -a 



1 d-1 
< — + 



d 

Y — 

where X^ and X$ are the arithmetic and harmonic means of a and f3 with 
proportions 9, 

X e =(- + ^--^) 1 and A+ = 9a + (l-9)/3. 



.a (3 

Proof. We will prove that Gg = Gg 
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We start proving that Gg C Gg and we consider A* € Gg. Therefore 
there exists a sequence of matrices {-A n }neiv of the form axn^d + /3(l — Xn)Id 
with Xn — v such that for every right-hand side /(x) the solutions u n of 

div (A n (x)V«„(a;)) = f(x) mD 
u n = on <9L> 

satisfy 

u n ^u w.Hl(D), 



A n Vu n -± A*Vu in L 2 (D) d , 
where u is the solution of the homogenized equation 

div (A*(x)Vu(x)) = f(x) in D 

u = on 9D. 

It is enough to observe that when f(x) is replaced by f(x,u) the above 
convergence holds P-a.e. w£0, and therefore A* £ Gg. 

We prove now that Gg C and we take A* € Gg. Therefore there 
exists a sequence of matrices {A n } n£ N of the form ax n I + — XnK with 
Xn — v such that for every right-hand side f(x,u) the solutions u n of 

div (4(i)Vii n (i,w)) = /(ic,o;) in D 

u n = on 3D { ' 

satisfy 

u n (-,oj) — in Hq(D), P-a.e. lo € . . 

AnVitnC-jW) A»Vu(-,w) inL 2 (D) a! , P-a.e. w € ' 

where u(-,uj) is the solution of the homogenized equation P-a.e. w E 

■div (A*(a;)V'5(a;, a;)) = f(x,oj) mD , , 

u = on <9L>. U 

Integrating the above expressions ([6]), 0, §E§ with respect to the random 
variable well and setting 

u(x) = / u(x,oo) dP(oj), u n (x) = I u n (x,oj) dP(u), 
Jn Jn 



and f(x) = / /(a;, w) dP{u), one has 



div(A„(x)Vu n (x)) = f{x) mD 

u n = on dD. 



and 

u n ^u in Hq(D), 



A n Vu n ->> A*Vu in L 2 (L>) d 
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where u is the solution of the homogenized equation 

-div (A*(x)S/u(x)) = f(x) in D 

u = on dD. 

/.From the generality of / we can deduce that A* G Gq. □ 
3.2 Relaxation 

We now consider the classical optimal design problem 



(O c ) minl(x) = / / j (x , u , u x (x , u)) dx dP(u) 



subject to 

X G {0, 1}), with A = aI dX + /3I d (l - X ), 

— div (j4(x)V«(x, oj)) = f(x,oj) in D, 
u = on <9D, 

P-a.e. w G f2, and to the volume constraint 

X(x) da; < L, 



D 

with L G (0, \D\). 

The lack of optimal solutions for problems of the type (O c ) is well known 
even in the deterministic case (see [12]). 

The basic idea for the relaxation process consists in considering a larger 
class of admissible designs, in order a new (relaxed) problem on this larger 
class admits optimal solutions. Having in mind the above Theorem [H we 
consider the space of generalized designs 

92) = {(#, ^4*) £ L°°(D, [0, 1]) x M dxd : A t eG %) a.e. x G Z?|. 

Therefore we define the relaxed version (O r ) of the above optimal design 
problem as 

(O r ) mm I (9, A*) = / / j(x, oj, u(x, oj)) dx dP(uj) 
subject to 

9 G L°°(D; [0, 1]), with A* G G e , 
-div (A*(x)Vu(x,u}j) = f(x,oj) in D, (9) 
u = on <9Z), 

P-a.e. w G fi, and the volume constraint 

9{x) dx < L. 

D 



7 



Theorem 2. (O r ) is a relaxation of (O c ) in the sense that 

1. the infima of both problems coincide 

2. there are optimal solutions for the relaxed problem (O r ). 

Proof. See for instance pQ Section 3.2, Theorem 3.2.1. □ 

4 Optimal solutions 

In this section we consider the above problems (O c ) and (O r ) in the special 
case when the cost functionals are either the compliance 

j(x,u,u) = f(x,u)u 

or the energy 

j(x,u,u) = -f(x,u)u 
and we will prove that in these situations our original design problem 



(O) min/(a) = / / j(x,ui,u a (x,u)))dx dP 
Jn Jd 

subject to, 

a € L°°(D), with a < a < j3, 
— div (a(x)Vu(x)) = / in D, 

u = on dD, 

f D a(x) dx = in 

admits optimal solutions. 

Theorem 3. In the cases either of the compliance or of the energy the 
optimization problem @ admits a solution. 

We know that the (O r ) problem admits optimal solution (since is the 
relaxed problem of (O c )), we check that these optimal solutions are solutions 
of our problem (O) from which we deduce the well-posed character of our 
problem. 

We analyze the optimality condition for (O r ) for the matrix A*. We 
denote by p the adjoint state, which is the unique solution in Hq(D) P-a.e. 
ui € of the adjoint state equation 

-div (A*(x)Vp(x,u))) = mD i (io) 

p = on dD, 

We remark that from (|10p it follows that in the compliance case we have 
p = u and in the energy case we have p = —u. 
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We fix the density G L°°(D; [0, 1]) and we introduce the Lagrangian 

L(M,(j),ip) = / / j(x,u,(j)(x,u}))dxdP(u)) + 

/ / [div (M(x)V</>(x,u;)) + f(x,cj)]ipdxdP(w) 
Jq Jd l - 1 

for any M G L°°(D;M dxd ) and 0,i/> G flJ(D) P-a.e. u G 0. We compute 
the partial derivative for any variable. 

It is clear that (J^L(M, </>, ip); = taking <fi = u solution of Q. We 

then determine the solution p so that, for all 4>\ G ^(D) P-a.e. w 6 H, we 
have 

—L(M,(j),p);<f)i) = 0, 



which leads to the formulation of the adjoint problem (|10p . 
Finally, from I(-,M) = L(M,u,p) it is easy to compute 

(dF £(M ' M ' p);Ml ) = ~ / / M iV«-VpcfedP(u;). 

Therefore it is easy to deduce that our cost functional is Gateaux differen- 
tiable with respect to the matrix variable and its derivative in the direction 
Mi is given by the above formula with u and p solutions of the state and 
adjoint equation respectively. Hence if M* if optimal, the optimality condi- 
tion 

() £(M*,u,p);M- M*\ > VM G G e 



,<9M 
becomes 

-/ f (M — M*) Vu • Vp dP(w) dx >0 VM G G$. (11) 
Jd Jn 

^From the optimality condition (|11|) we obtain 
/ / M*Vu- VpdP(u))dx = max / / MVu ■ VpdP(w) dx. (12) 

J D J IS/1 £l(jt Q J J 

Using the algebraic expression 

AMy ■z = M{z + y)-{z + y)- M(z - y) ■ {z - y), 
we have that for any M G Gg 

4MVu • Vp < max A(Vu + Vp) ■ (Vu + Vp) - min A(Vu - Vp) ■ (Vu - Vp) 

and using the identification of Ge we obtain 

4MVu • Vp < At|Vu + Vp| 2 - A„ |Vu - Vp| 2 . 
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Having in mind that in our problem u = p for the compliance and u = —p 
for the energy, the necessary condition above reads 

/ / M*V«- VudP{uj)dx = / / Xg\Vu\ 2 dP{uj) dx (13) 
Jd Ja J d Jn 

with Xg = A^~ for the compliance and Xg = Xg for the energy. 

If we analyze the optimality condition (|13p we obtain that Xg is an eigen- 
value of M* and Vu(x,u) is an eigenvector P-a.e u S Q, i.e., 

M*Vu = XgVu P-a.e. we!] (14) 

where u is the solution of the state equation. 
For the compliance case u = p and 

M*Vu • Vu = A^|Vn| 2 , 

there exist several matrices in Gg with this property, and it is enough to 
take a rank one laminate with normal direction of lamination n orthogonal 
to Vu and the optimal volume fraction A^" . 
For the energy case u = —p and 

M*Vu ■ Vu = A^|Vn| 2 , 

there exist an unique matrix in Gg with this property which corresponds to 
a rank one laminate with normal direction of lamination n parallel to Vit 
and the optimal volume fraction A^. 

We would like to stress that for any case the optimal matrix M* is a first 
order laminate with deterministic optimal volume fraction Xg and random 
direction of lamination according with the random value of Vu. 

Therefore, from the analysis of the optimality condition we concluded 
that the optimal matrix M* € Gg verifies the condition (|14p . We remark that 
this condition does not imply that the optimal matrix M* is the isotropic 
matrix Xgld', the important conclusion of (j!4[) is that Xg is an eigenvalue of 
M* and Vu is an associated eigenvector, where u is the solution of the state 
equation. In particular, this implies that the optimal value from the (O r ) is 
attained on the simpler problem 

(O) min/(#) = / / j(x,cj,u(x,u)))dx dP(u) 
JuJd 

subject to, 

9eL°°(D; [0,1]), 
- div (Xg(x)Vu(x)) = f inD, 
u = uo on dD, 
f n 9{x) dx = L 

Finally, the proof of Theorem reduces to the fact that for 6 6 L°°(D; [0, 1]) 
one has that Xg S L°°{D; [a, /3]), in particular it is enough to take a = Xg to 
deduce the existence of optimal solution of our original problem (O). 
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5 Numerical analysis of the optimal design prob- 
lem 



We approach in this section the numerical resolution of the problem (O) 
for the compliance case for which j(x,uj,u) = u(x)f(x,uj) with Dirichlet 
boundary condition, for which the cost can be written as 

a(x)\Vu(x,u))\ 2 dxdP(uj). 



n jd 

We first describe an algorithm of minimization and then present some nu- 
merical experiments. We treat the problem 

(O) min/(a)= f [ u(x)(f(x) + Z(x,u))dP(u)dx 
Jd Jn 

subject to, 

a € L°°(D), with a < a < /3, 

(15) 



a € L°°(D), with a < a < /3, 
- div (a(x)Vu) = / + £ in D, P-a.e. u6i! 
■u = on dD, 
J D a(x) dx = m 

where ^ is a random variable £ = £(x,w). Similar computations are also 
made for the energy case j{x,uj,u) = —u(x)f(x,uj). 

5.1 Algorithm of minimization 

We present the resolution of the optimal design problem (O) using a gradient 
descent method. In this respect, we compute the first variation of the cost 
function with respect to a. 

For any r/ € M + , rj -C 1, and any a € L°°(D), we associate to the 
perturbation aJ 1 = a + rjd of a the derivative of 7 with respect to a in the 
direction a as follows: 

07(a) , 7(a + ria) - 7(a) 
— - — • a = hm . 

da r?— >o 7] 

Theorem 4. The first derivative of I with respect to a in any direction a, 
exists and takes the form 



07(a) 
da 



a = ~J ®(J VuVpdP(w)) dx, (16) 



where u is the solution of (|15p and p is the solution in Hq(D) of the adjoint 
problem 

— div (a(x)Vp) = / + £ in 7>, P-a.e. uj € , . 

p = on dD. 
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Proof. We introduce the Lagrangian 



L(a,4>,ifj) = / / 4>(x, ui)(f(x) + £(x, uj)) dP(u) dx + 
J d Jn 

div (a(x) V(j)(x, u)) + /(x) + ^ (x, cj) i/> dx dP(u) 




d Jn 



for any a G L°°(D; [a,0\) and <fi,ip E Hq{D) P-a.e. We compute the 

partial derivative for any variable. 

It is clear that {-§^£>{a, (/), tp); ^i) =0 taking cj) = u solution of (fT5j) . We 
then determine the solution p so that, for all fa G H l (D), we have 







L(a,(j),p);(f)i\ = 0, 



which leads to the formulation of the adjoint problem (|17p . 
Finally, it is easy to compute 

(J^L{a,u,p)-,aJ = - J a(J VuVpdP(u;)) dx. (18) 

Next, writing that 1(a) = £>(a,u,p), we obtain ()16p from (|18j), □ 

In order to take into account the volume constraint on a, we introduce 
the Lagrange multiplier 7 G R and the functional 

7 7 (a) = /(a) + 7 / a(x)dx. 



'7 

Using Theorem HI we then obtain easily that the first derivative of I~ is 
d/ 7 (a) _ _ _y y u y pdp ^ jJ ^ dx + 1 J a(x)dx, 

which leads us to define the following descent direction: 

a(x) = (y VuVpdP(u) - 7) Vx G £>. (19) 

In this way, for any function 7/ G L°°(D,M + ) with || 77 1| 2,00 m) small enough, 
we have I 7 (a + 77a) < I 7 (a). The multiplier 7 is then determined so that, 
for any function r\ G L°°(Z),R + ), ||a + 7?a||Li(D) = m ! leading to 



(20) 






f a(x)dx — rn^j 


+ J 




' VuVpdP(uj)dx 
n 


J 


/ rj(x) dx 

D 
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where the function rj is chosen such that a + r/a G [a, /3] for all x G D. A 
simple and efficient choice consists of taking rj(x) = e(a(x) — a)(f3 — a(x)) 
for all x € D with e small and positive. 

Finally we show the descent algorithm to solve numerically the optimiza- 
tion problem (O). 

We consider < a < (3, m € P\D\) and e < 1, E\ -C 1 data of the 

problem, the structure of the algorithm is as follows. 

• Initialization of the density a S L°°(D; [a,/3]); 

• for k > 0, iteration until convergence (i.e., \Iy(a k+1 ) — I^{a k )\ < 
£i|/ 7 (a°)|) as follows: 

— compute the solution u a k of (fT5|) and then the solution p a k of 
(fT7|) . both corresponding to a = a k ; 

— compute the descent direction a defined by (fT9j) . where the mul- 
tiplier 7 is defined by (p0|) : 

— update the density a k in D: 

a k+1 = a k + e(a k -a){P -a k )a k , 

with e € M + small enough to ensure the decrease of the cost 
function, a k+1 G L°°(D, [a, (3}). 

5.2 Numerical experiments 

In this section we implement the gradient descent algorithm explained in 
the previous subsection. These sort of problems have been studied by other 
authors ([I], [3], [5], [IS]), and we will show the numerical results according 
with previous numerical simulations of the previous authors. 

We solve the problem (O) on the square domain D = (0, l) 2 for two 
phases a = 1 and (3 = 2, the determinist part of the right-hand side of the 
state equation / = 1 and we consider the volume constraint m = = 1.5, 
i.e. we can use the same amount of a or /3 mass. In order to simplify the 
numerical computations we choose the random variable £ with a discrete 
distribution of probability. We consider two different cases for £: 

• Case 1: = ± X d where D = [±, | ] 2 C D 

• Case 2: = ±x.Di where D\ = D \ Dq 

and in both cases -P({£ = x}) = = ~x}) = I- It could be possible 

that the algorithm fall down into local minima of /, for this reason, we 
consider a constant initialization ao = m, without any privilege for the 
optimal localization. 
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In the figures below we represent the corresponding optimal mass distri- 
bution a. We show numerical results for full deterministic case and the two 
random cases described above, all for the compliance and energy minimiza- 
tion. The result are qualitatively different for any case. 



Optimal design determinist Optimal design determinist 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 1: Optimal distribution for the full deterministic case. Left: Com- 
pliance minimization. Right: Energy minimization. 



Optimal design random Optimal design random 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 2: Optimal distribution for the random Case 1. Left: Compliance 
minimization. Right: Energy minimization 

With respect to the compliance minimization, we observe that the limit 
densities follow a similar distribution, where the smaller amount of mass is 
placed as a cross-shape. Taking as reference the picture of Figure Q] (the 
deterministic case) we observe different densities for the two random cases. 
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Optimal design random 



Optimal design random 



0.3 



0.1 



0.2 



0.5 



0.6 



0.9 



0.8 



0.7 



0.4 




0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



Figure 3: Optimal distribution for the random Case 2. Left: Compliance 
minimization. Right: Energy minimization 

For the Case 1, where the random perturbation is at the middle of the square 
a bigger amount of mass (the white color at the pictures) is distributed for 
this place and the thickness of the cross is bigger at the corners (in order 
to save the volume constraint) (see Figure [2]). For the Case 2 the effect is 
the inverse, in this case the random perturbation is near the boundary; the 
optimal distribution consists in a smaller amount of mass in the middle and 
a bigger concentration near to the boundary with a smaller thickness of the 
cross (see Figure ED . 

For the energy minimization, the optimal distribution is fully different 
and it is in accord with previous analysis (see [I], [5]). We take again the 
deterministic case as the reference design; for this case the simulations show 
that a bigger amount of mass is placed at the corners and at a square in 
the middle of the domain of design. For the Case 1 the optimal distribution 
is very similar to the deterministic one, and the changes correspond with 
placing a bit more of mass at the corners of the domain. For the Case 2, 
the effect is the reverse, in this case there are less mass at the corners and 
almost all the mass is placed at a square in the middle of the domain. 

Finally, and in conclusion the numerical experiments indicate that under 
random forces the optimal distribution consists in placing a bigger amount 
of mass where this random force acts in, in order to take into account the 
possibility for the load to have a random variation. 
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